Surrogate model for an aligned-spin effective one body waveform model of binary neutron star inspirals using Gaussian process regression

Fast and accurate waveform models are necessary for measuring the properties of inspiraling binary neutron star systems such as GW170817. We present a frequency-domain surrogate version of the aligned-spin binary neutron star waveform model using the effective one body formalism known as SEOBNRv4T. This model includes the quadrupolar and octopolar adiabatic and dynamical tides. The version presented here is improved by the inclusion of the spin-induced quadrupole moment effect, and completed by a prescription for tapering the end of the waveform to qualitatively reproduce numerical relativity simulations. The resulting model has 14 intrinsic parameters. We reduce its dimensionality by using universal relations that approximate all matter effects in terms of the leading quadrupolar tidal parameters. The implementation of the time-domain model can take up to an hour to evaluate using a starting frequency of 20Hz, and this is too slow for many parameter estimation codes that require $O(10^7)$ sequential waveform evaluations. We therefore construct a fast and faithful frequency-domain surrogate of this model using Gaussian process regression. The resulting surrogate has a maximum mismatch of $4.5\times 10^{-4}$ for the Advanced LIGO detector, and requires 0.13s to evaluate for a waveform with a starting frequency of 20Hz. Finally, we perform an end-to-end test of the surrogate with a set of parameter estimation runs, and find that the surrogate accurately recovers the parameters of injected waveforms.Comment: 19 pages, 10 figures, submitted to PR

Aligned spin neutron star-black hole mergers: a gravitational waveform amplitude model

The gravitational radiation emitted during the merger of a black hole with a neutron star is rather similar to the radiation from the merger of two black holes when the neutron star is not tidally disrupted. When tidal disruption occurs, gravitational waveforms can be broadly classified in two groups, depending on the spatial extent of the disrupted material. Extending previous work by some of us, here we present a phenomenological model for the gravitational waveform amplitude in the frequency domain encompassing the three possible outcomes of the merger: no tidal disruption, "mild" and "strong" tidal disruption. The model is calibrated to 134 general-relativistic numerical simulations of binaries where the black hole spin is either aligned or antialigned with the orbital angular momentum. All simulations were produced using the SACRA code and piecewise polytropic neutron star equations of state. The present model can be used to determine when black-hole binary waveforms are sufficient for gravitational-wave detection, to extract information on the equation of state from future gravitational-wave observations, to obtain more accurate estimates of black hole-neutron star merger event rates, and to determine the conditions under which these systems are plausible candidates as central engines of gamma-ray bursts, macronovae and kilonovae.Comment: 15 pages, 7 figures, 1 tabl

Effective-one-body waveforms for binary neutron stars using surrogate models

Gravitational-wave observations of binary neutron star systems can provide information about the masses, spins, and structure of neutron stars. However, this requires accurate and computationally efficient waveform models that take <1s to evaluate for use in Bayesian parameter estimation codes that perform 10^7 - 10^8 waveform evaluations. We present a surrogate model of a nonspinning effective-one-body waveform model with l = 2, 3, and 4 tidal multipole moments that reproduces waveforms of binary neutron star numerical simulations up to merger. The surrogate is built from compact sets of effective-one-body waveform amplitude and phase data that each form a reduced basis. We find that 12 amplitude and 7 phase basis elements are sufficient to reconstruct any binary neutron star waveform with a starting frequency of 10Hz. The surrogate has maximum errors of 3.8% in amplitude (0.04% excluding the last 100M before merger) and 0.043 radians in phase. The version implemented in the LIGO Algorithm Library takes ~0.07s to evaluate for a starting frequency of 30Hz and ~0.8s for a starting frequency of 10Hz, resulting in a speed-up factor of ~10^3 - 10^4 relative to the original Matlab code. This allows parameter estimation codes to run in days to weeks rather than years, and we demonstrate this with a Nested Sampling run that recovers the masses and tidal parameters of a simulated binary neutron star system.Comment: 17 pages, 11 figures, submitted to PR

Alignment and Aperture Scan at the Fermilab Booster

The Fermilab booster has an intensity upgrade plan called the Proton Improvement plan (PIP). The flux throughput goal is 2E17 protons/hour, which is almost double the current operation at 1.1E17 protons/hour. The beam loss in the machine is going to be the source of issues. The booster accelerates beam from 400 MeV to 8 GeV and extracts to the Main Injector. Several percent of the beam is lost within 3 msec after the injection. The aperture at injection energy was measured and compared with the survey data. The magnets are going to be realigned in March 2012 in order to increase the aperture. The beam studies, analysis of the scan and alignment data, and the result of the magnet moves will be discussed in this paper.Comment: 3 pp. 3rd International Particle Accelerator Conference (IPAC 2012) 20-25 May 2012, New Orleans, Louisian

Mars mission solar array Semiannual progress report, period ending 31 Dec. 1969

Design and testing of beryllium-structure solar panel for Mars missio